Self-organized criticality and interface depinning transitions 
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We discuss the relation between self-organized criticality and depinning transitions by mapping 
sandpile models to equations that describe driven interfaces in random media. This equivalence 
yields a continuum description and gives insight about various ways of reaching the depinning 
critical point: slow drive (self-organized criticality), fixed density simulations, tuning the interface 
velocity (extremal drive criticality), or tuning the driving force. We obtain a scaling relation for the 
correlation length exponent for sandpiles. 
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I. INTRODUCTION 

Many systems respond to external perturbations by 
avalanches which behave intermittently with a power-law 
distribution of sizes. The paradigm of such self-organized 
critical (SOC) behavior is the so-called sandpile model 
E|. It maintains by an infinitely slow drive a critical 
steady-state, where the internal dissipation balances the 
external drive. Candidates for such phenomena include 
granular piles [§J|, microfracturing processes Q, and 
earthquakes [|| . Despite many theoretical and numerical 
investigations a thorough understanding of self-organized 
criticality is still lacking ]6j-|l3|]. Fundamental problems 
which need to be solved involve deriving a continuum the- 
ory which would for instance determine the upper critical 
dimension, above which mean-field theory applies PJTif. 

Similar behavior can be found from elastic interfaces 
driven through random media JT5|-p^|. They undergo a 
continuous (critical) depinning transition as the exter- 
nal driving force is varied. With increasing force one 
passes from a phase where the interface is pinned to a 
depinned phase where the interface moves with a con- 
stant velocity. Close to the critical point, the motion of 
the interface takes place in "bursts" with no character- 
istic size and the interface develops scaling described by 
critical exponents. These phenomena can be met in flu- 
ids driven through porous media Jl9] , in domain walls in 
magnets (the Barkhausen effect) pOj , in flux lines in type 
II superconductors plj] , and in charge-density waves [B2| . 

In this paper we investigate the connections be- 
tween self-organized criticality and depinning transitions 
fl[R|il j20|Jp f6|. We first establish a generic, exact 
relation J26{] between sandpile models and driven inter- 
faces which builds upon previous investigations of e.g. a 
charge-density wave model Q and a rice-pile model Q . 
Specifically, we discuss the Bak-Tang-Wiesenfeld (BTW) 
[ j| model and, as an example, a stochastic sandpile model 
[ 2^,^8| through a mapping to a model for interface depin- 
ning with slightly different noise terms. 

The mapping enables one to understand the slow-drive 



criticality used in sandpile simulations in terms of stan- 
dard concepts for driven interfaces. Using the continuum 
theory for interface depinning it follows for these sand- 
pile models that the upper critical dimension d c is 4, and 
the relevant noise is of quenched type. The connection 
with interfaces allows us to establish a scaling relation for 
the correlation length exponent for sandpile models. In 
addition, we discuss in the interface representation sand- 
piles driven at fixed density, driven at boundaries, and 
extremal drive criticality. 



II. SANDPILES 

The sandpile models are here defined as follows: to 
each site of a d-dimensional lattice (square in d — 2) of 
size L d is associated a variable z x which counts the num- 
ber of grains on that site. When the number of grains 
on a site exceeds a critical threshold z x > z c , the site 
is active and it topples. This means that 2d grains are 
removed from that site and given to the 2d nearest neigh- 
bors (nn): z x — > z x — 2d, z nn — > z nn + 1, Vnn. Sandpiles 
are usually open such that grains which topple out of the 
system are lost (in one dimension: zq = zl + \ = 0). It is 
also possible, as discussed later, to use periodic bound- 
ary conditions. When there are no more active sites in 
the system, one grain is added to a randomly chosen site, 
z x ~ * %x + 1- The time and number of topplings till the 
system again contains no active sites define an avalanche 
and its internal lifetime and size. For the BTW model 
one has Zr = 2d — 1, whereas for stochastic sandpile 
models J27],|28| the threshold z c is not constant. Below 
we will focus on 1) the BTW model and 2) a stochastic 
model where the threshold z c is randomly chosen to be 
for example 2c? — 1 or 2c? after each toppling, i.e. P(z c ), 
the probability distribution of the z c 's, is any reasonable 
choice (i.e. decaying sufficiently fast). 

In terms of the internal avalanche time, the external 
drive is infinitely slow |2{|. After a transient, the system 
reaches a steady-state in which the slow drive and the 
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dissipation of grains balance each other. The boundary 
conditions (BCs) are essential to obtain criticality and 
they are usually of the Dirichlet type,z = 0, such that 
particles are dissipated to the outside . Alternatively, 
the SOC steady state can be reached by using bulk dis- 
sipation and, e.g., periodic BCs ylj. 

In the SOC steady-state the probability to have 
avalanches of lifetime t and size s follow power-law dis- 
tributions: p(t) = t- Tt f t (t/L z ) andp(s) = s- T f(s/L D ), 
with s ~ t D / z and z(r t - 1) = D(t-1) pj. Here the size 



1 /v = 2. 



(3) 



scales as 



and the (spatial) area as I (for compact 



avalanches) with I the linear dimension. The fact that 
each added grain will perform of the order of L 2 topplings 
before leaving the system leads to the fundamental result 



(s) ~ L 2 



(1) 



independent of dimension J^H. Thus, r = 2 — 2/D and 
r t = 1 + (D - 2)j z. Equation (0) yields 7/1/ = 2, 
where 7 describes the divergence of the susceptibility 
(bulk response to a bulk field) near a critical point, 
X = (s) ~ |A|~ 7 , and v is the (spatial) correlation length 
exponent, £ ~ |A| _iy Q. Here A = ( — ( c is the control 
parameter, ( — (z x ), and the critical value £ c = (z x ) soc , 
where this average is taken in the slowly driven SOC 
steady-state with A = 

III. INTERFACE DEPINNING 

For driven interfaces in random media critical scaling 
is obtained with a force F close to a critical value F c . 
Depinned interfaces move with a velocity v ~ f e , with 
/ = F — F c > 0. Pinned interfaces are blocked by pin- 
ning paths/manifolds which arise from the quenched dis- 
order environment. Close to criticality, correlations scale 
as x 2x , with x the roughness exponent, up to a corre- 
lation length £ ~ |/| _!/ - The characteristic time scale 
is £ z with z the dynamic exponent and it follows that 
9 = v{z — x) |l5| , |l6| ]. Near the depinning transition, the 
simplest choice to describe the dynamics of the interface 
is the following continuum equation ('quenched Edwards- 
Wilkinson', or linear interface model, LIM) fl^ |l7||: 



dH 
~dt 



V 2 ff + r)(x,H) +F. 



(2) 



Here, H (x, t) measures the height of a given site x at time 
t. The quenched noise 77(0;, H) has correlations given 
by (r)(x,H)r](x',H')) = S d (x - x')G{H - H'), where 
G(H — H') decays rapidly, approximated by a delta func- 
tion for random- field disorder. The critical exponents 
at the depinning transition have been calculated by e- 
expansions |l5|,[l6| and simulations [ p7p^j3l[ |. The up- 
per critical dimension is d c — 4, above which mean-field 
theory applies p2| . Below we will also discuss so-called 
columnar noise with G(H) = 1 |33|] . 

The interface equation (^|) obeys an invariant so that 
the static response scales as xil^ = 0) ~ q~ 2 , i.e., 



For forces below F Cl the (bulk) response of the inter- 
face triggered by a small increase in F scales as xtmik = 
d {H)/dF - (F c - F)~ 7 . Right at the critical point one 
can argue as follows |L5|,|l6|: the roughness of the inter- 
face scales as l x and assuming that A (H) will scale in 
the same way it follows 



7 = 1 + X v - 



(4) 



This yields x + V^ = 2, i.e., there are only two inde- 
pendent exponents for depinning described by (^|). The 
standard scaling relations are valid for interfaces with 
parallel dynamics: all sites with dH/dt > are updated 
in parallel. Note that interfaces with extremal (i.e., one 
unstable site at a time) and parallel drive have the same 
pinning paths. This manifests the Abelian character of 
the LIM in that the order in which active sites are ad- 
vanced does not matter fE4|. 



IV. MAPPING OF SANDPILE DYNAMICS 

Next we will show that the SOC critical behavior can 
be related exactly to the slowly driven depinning transi- 
tion in an interface model. Thus, Eqs. (lj) and (||) are 
equivalent and Eq. (|i|) yields an expression for the cor- 
relation length exponent v for sandpiles. The first step 
is to formulate the stopping of an avalanche in a SOC 
system as being due to a pinning path for an interface 
H(x,t). This field is given in the continuum limit by 



H(x,t) 



dt' p(x,t'), 



(5) 



where the order parameter p(x, t) is the activity (top- 
plings) at site x at time t, i.e., p = H = v ~ f e . In 
words: H {x 1 1) counts the number of topplings at site x 
up to time t. At the microscopic level this is an exact 
correspondence between a toppling and the interface ad- 
vance. A toppling takes place when z x > z Cl which by 
the relation 



dH 



(6) 



yields the dynamics dH/dt >0 => H — > H + l, whereas 
H is unchanged at the sites where no toppling takes place. 
The dynamics of sandpile models thus map to discrete 
interface equations where an avalanche takes the inter- 
face H (x, t) from one pinning path to the next in the 
quenched random medium |l^,|l^,^,^3|,^5| . Since the in- 
terface counts topplings it does not move backwards and 
thus Eq. (^|) effectively reads dH/dt — 8(z x — z c ), which 
is the standard discretization for depinning models p7| . 
We are currently investigating the applicability of such 
discretization procedures to various models pqBJ]. 



2 



Next, we express z x in terms of H (x, t) for the specific 
models introduced above. The number of grains z x on 
site x is z x = Ni n — N out + F(x,t), where iVj„ is the num- 
ber of grains added to this site from its 2d nearest neigh- 
bors (nn) and N out is the number of grains removed from 
this site due to topplings. The (external) driving force 
F (x, t) counts the number of grains added from the out- 
side. Since N in = E nn H(x nn ,t) and N out = 2dH{x,t) 
(for details and extensions to other models see p4|) we 
arrive at 



dH 
~dt 



V 2 H - z c (x,H) + F(x,t), 



(J) 



where V 2 H is the discrete Laplacian. The Dirichlct 
boundary conditions for z x become H = and the dy- 
namics is parallel. Similar connections have been previ- 
ously discussed for a charge-density wave model |23] and 
for a boundary driven rice-pile model flfl] (see below). 
In the stochastic model, z c (x,H) is a random variable 
which changes after each toppling. Thus z c (x, H) acts 
like quenched random point-disorder similar to r](x, H) in 
Eq. (^). The BTW model has z c equal to a constant. The 
dissipation needed to reach the SOC state (loss of grains 
z x ) takes place through the BC of H = 0. Using strong 
boundary pinning may thus give rise to the possibility 
of observing SOC experimentally in systems displaying 
a depinning transition. Wc emphasize that the mapping 
prescription can in principle be applied to any sandpile 
model. For other, more complicated, toppling rules jp|]35| 
additional terms like the "Kardar-Parisi-Zhang" nonlin- 
earity |V-ff| 2 may appear. 

On the internal (fast) time scale the driving force 
F(x, t) does not act as a time-dependent noise but as 
columnar- type disorder. It counts all the grains added to 
the system by the slow drive, i.e. F(x,t) — » F(x,t) + 1, 
and thus increases as function of time in an uncorrelated 
fashion. In the opposite limit when a grain is added (e.g.) 
each time step ("fast drive") F(x, t) would correspond to 
a time-dependent noise Since H = at the bound- 
ary and F increases as function of time the steady-state 
profile of H will be close to a paraboloid or, in one dimen- 
sion, a parabola (see also pH). In the steady-state, just 
after an avalanche, the slowly increasing force F is bal- 
anced by the negative curvature ~S7 2 H of the paraboloid 
such that all sites are pinned (dH/dt < 0). This illus- 
trates that the interface effectively is driven by a force 
equal to the critical force F c = Q c — ~z~ c , where ~ c is the 
average of z c (x, H) in the steady state (for the BTW 
model trivially z c = 2d — 1). Accordingly, the slow drive 
reaches the depinning critical point by adjusting the dis- 
sipation to the driving force such that the velocity (order 
parameter) is infinitesimal. 

The steady-state of the different sandpile models is de- 
scribed by an equation similar to Eq. (g). Thus the ex- 
ponent relation (||) holds and it is equivalent to Eq. (|l|) 
which describes the scaling of the average avalanche size 
( "susceptibility" ) . Assuming that a roughness exponent 
X can be defined for sandpile models, one can argue that 



Eq. (0) is valid also for sandpiles. Furthermore, the upper 
critical dimension is d c — 4. Note that the ensuing noise 
will contain a columnar component |2^,^5|,[33| due to the 
random drive F(x,t). The one-dimensional BTW model 
has a critical force F c = 1 — 1 = 0, which corresponds 
to the critical point of the columnar-disorder interface 
model j33|. In d > 1, one has F c < which in combina- 
tion with the fact that the interface by definition cannot 
move backwards implies that the BTW model displays 
a more complicated behavior than the columnar mod- 
els investigated in [^3 33 1. Note also that avalanches in 
stochastic models will have a random structure due to the 
explicit point disorder whereas avalanches in the BTW 
model show a more regular behavior p6| , [36| |. 

For the case of the boundary driven one-dimensional 
rice-pile models [ p7[p8| a similar mapping of the dynam- 
ics can be done with an auxiliary field H(0,t) and a 
drive implemented as H(0,t) -> H(0,t) + 1 Q. The 
rice-pile models have Dirichlet BC at x = and Neu- 
mann BC (reflective) at x = L which yields (s) ~ L. In 
our picture the boundary drive is F(l,t) — ► F(l,t) + 1 
and F(x > l,t) = 0. Because of the Neumann BC 
\H(L,t) = H(L + l,t)] the steady state develops a 
parabolic profile with the left branch pointing up [fLOf. 



V. VARIOUS ENSEMBLES 



We next consider the more straightforward cases in 
which sandpiles are studied with periodic boundary con- 
ditions (amounting to H(l) = H(L) in one dimension). 
In such cases the SOC steady state can be tuned into by 
various approaches. It can be reached by using a care- 
fully tuned bulk dissipation e ~ L~ 2 |13[ . In this case, 
periodic BCs are also the best since the scaling of the 
system is not a mixture of boundary and bulk scaling 
|p7| . As above, we arrive at 

BJ-f 

— =\7 2 H-z c (x,H)-e(x,H)+F(x,t) (8) 
at 

with H periodic. As in Eq. ([?]), the force F(x, t) is colum- 
nar and increases on the slow time scale. The dissipation 
e(x, H) takes now into account all the grains removed 
before the site at x topples. It increases with a (small) 
probability only when a site topples and this means that 
e explicitly depends on H . Therefore, a dissipation event 
effectively corresponds to a shift in the z c value. Thus, 
one obtains that the BTW model with bulk dissipation 
contains a very weak point-disorder component (since the 
increases in F equal in the statistical sense the increases 
in e). Though point-disorder is in general expected to be 
a relevant perturbation, in the infinite system size limit 
the Larkin length Jl^,[l6| associated with the cross-over 
from columnar behavior diverges and thus the avalanche 
behavior is not governed by the weak point disorder. By 
this argument the BTW models with or without bulk 
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dissipation are equivalent to the same interface depin- 
ning equation (|^) in accordance with simulations of the 
BTW and bulk dissipation models |Q. Note that the 
boundary critical behavior of the BTW model depends 
on the specific boundary condition: Dirichlct BCs display 
a different behavior 39 1, whereas Neumann BCs (reflec- 
tive) are similar to the bulk. In the case of periodic BCs 
and bulk dissipation, the ii-field fluctuates around an 
average flat profile. The terms F(x,t) and e(x,H) will 
balance each other in the steady state with an average 
difference such that F c — £ c — < 0. For larger dissipa- 
tion rates the system moves away from the critical point 
and, in analogy to Eq. (^), the bulk susceptibility scales 
as Xbuik ~ l/e ~ with v t = 1/2 [p[. 

The fixed density (or energy) drive previously used in 
simulations |j 11 corresponds to a normal driven inter- 
face. Thus, the situation is such that H(x,t = 0) = 0, 
C = L- d J2 X F(x, 0) with F(x, t) = F(x, 0), and periodic 
BCs and e(x,H) = such that no 'grains' are lost. The 
control parameter A = £ — £ c (= F — F c = /) is varied 
and criticality is only obtained when A = 0; note that 
choosing £ corresponds to using a spatially dependent 
force F(x,0) with ( = (F(x,0)). Here, the system is 
not generally in the SOC steady-state but by letting the 
control parameter A — > one reaches the critical point 
H[ll] . The noise is set at the beginning of an avalanche 
at the columnar values F(x,0). Depending on the ex- 
act nature of the initial configuration one may observe a 
different dynamic behavior but the steady-state behavior 
should correspond to the slowly driven case jll| . 

In "microcanonical" simulations one has dissipa- 
tion operating on the slow time scale with exactly the 
same rate as F(x,t). Thus microcanonical simulations 
correspond to fixed density simulations with a specific 
initial configuration: after each avalanche, the time is re- 
set to zero, the force is replaced with F — ► F + V 2 H, and 
the forces at x' (x") are increased (decreased) by one unit 
where x' and x" are randomly chosen sites. Finally the 
interface is initialized, H = 0. Since the V 2 H term does 
not introduce correlations this new starting condition is 
equivalent to the fixed density case but with the initial 
configuration chosen to be in the SOC steady state. 

Combining the scaling relations (|3|) and (|J) it follows 
that 



2 + d= D + l/v, 



(9) 



where D = d + \. In addition, the average area scales 
as (l d ) ~ L x l v . These relations are also valid for sand- 
piles and Eq. ([)]) provides estimates for v. in d = 1, 
v w 1.30, and in d = 2, v w 0.78. Numerical results 
yield v = 1.25(5) (d = 1, stochastic model) |4l[] and 
v = 0.79(4) (d = 2, BTW model) 0. Note, however, 
that the estimates quoted for v for sandpile models de- 
pend on the relation D = d + x, which means that the 
underlying assumption is that the roughness exponent x 
can be defined for slowly driven sandpile models. 



VI. CONCLUSIONS 



In summary, we have started from the depinning equa- 
tion (||) to discuss the continuum description of self- 
organized critical sandpile models. Thus, their upper 
critical dimension is d c = 4 and a scaling relation for 
the correlation length exponent v is obtained. We find 
that the BTW model has columnar disorder F on the 
avalanche time scale whereas the stochastic models have 
explicitly point disorder included. Other models with 
slightly modified toppling rules (e.g., the Manna model 
p2| ) may or may not belong to the same classes depend- 
ing on the noise terms arising from the mapping (this 
we are currently investigating further in p6| , f34| ) . The 
present approach shows that the relevant noise for sand- 
piles is 'quenched'. The physics of sandpiles is such that 
the random decisions or events (grain deposition, choices 
for thresholds) are frozen into the dynamics of a site as 
long as it is stable, and their memory decays only slowly 
as the activity goes on. A recent field theory for p(x,t) 
used analogies from systems with absorbing states and 
assumed that the noise was Reggeon field-theory like (i.e., 
time-dependent and not quenched) Physically, the 

effect which is not incorporated in such Gaussian cor- 
relations is that the pinning forces along the interface 
selects a pinning path in the random media which stops 
the avalanche. 

The mapping between interface and sandpile dynam- 
ics allows one to characterize the sandpile universality 
classes by the quenched noise in the interface equations. 
It also allows to gain novel insight about the previously 
introduced ways of reaching the depinning critical point: 
balancing the force with dissipation (slow drive, or self- 
organized criticality), tuning the average force (as for 
fixed density sandpiles) , tuning the interface velocity (ex- 
tremal drive criticality), and finally tuning the driving 
force. This becomes possible because of the diffusive 
character of interface or sandpile dynamics and because 
of the Abelian character of the linear interface equation. 

K. B. L. is supported by the Carlsberg Foundation. 
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